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Abstract 

The shower simulation code CORSIKA has been used to investigate the electron en- 
ergy and angular distributions in high-energy showers. Based on the universality of 
both distributions, we develop an analytical description of Cherenkov light emission 
in extensive air showers, which provides the total number and angular distribution 
of photons. The parameterisation can be used e.g. to calculate the contribution 
of direct and scattered Cherenkov light to shower profiles measured with the air 
fluorescence technique. 
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1 Introduction 



The measurement of ultrahigh-energy air showers above 10 eV is one of the 
main tools to understand the nature, origin and propagation of cosmic rays 
at highest energies. Large-scale experiments, like the Pierre Auger Observa- 
tory [1], AG ASA [2], HiRes [3], and Telescope Array [4] focus on the precise 
determination of the energy spectrum, mass composition and arrival direction 
distribution of ultrahigh-energy cosmic rays (UHECR). 

The calorimetric measurement of the longitudinal shower profile with the flu- 
orescence technique is one of the most direct methods to determine the pri- 
mary energy [5] and also can be used to infer the primary particle type. The 
charged shower particles induce nitrogen fluorescence light, which is emitted 
isotropically mainly in the near UV-range of 300-400 nm [6]. Assuming the 
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fluorescence yield is proportional to the energy deposit by a shower, as indi- 
cated by measurements [7,8], observing the emitted fluorescence light gives a 
measure of the local ionisation energy deposit and hence the shower size. 

However, charged particles in extensive air showers (EAS) are mostly relativis- 
tic and Cherenkov light is produced as well. In contrast to the fluorescence 
light, Cherenkov photons are emitted mostly in the forward direction. Depend- 
ing on the observation angle with respect to the shower axis, the intensity of 
produced Cherenkov light contributes a non-negligible signal in fluorescence 
detectors [5,9,10,11], which can be of same order as the fluorescence light 
signal itself [10,11,12]. Thus, an efficient correction for or an explicit consider- 
ation of the Cherenkov contribution to the measured signal is needed for the 
determination of the primary particle properties. 

In the pioneering Fly's Eye experiment, Cherenkov light was estimated [5] 
based on Hillas' simulations of 100 GeV photon showers [13]. Assuming uni- 
versality of electron distributions in EAS, the parameterisation in shower age 
developed by Hillas was applied to showers of ultrahigh-energy [5,14,15]. More 
recent simulation studies at higher energies and with hadron primaries confirm 
such an universality of electron energy spectra [16,18]. 

During the last decades substantial progress has been made in the simulation 
of very high energy EAS. New simulation packages such as CORSIKA [19] 
and modern high-energy interaction models allow a more reliable prediction 
of particle distributions in showers. At the same time the precision of UHECR 
detectors has increased dramatically, requiring a treatment of Cherenkov light 
production as accurate as possible. 

In this work we perform a systematic study of shower particle distributions 
relevant to Cherenkov light calculation. A complete model for Cherenkov light 
calculation is developed and compared to predictions from CORSIKA employ- 
ing QGSJET01 [20]. In this approach, the total Cherenkov light can be calcu- 
lated as function either of the shower size profile or the energy deposit profile. 
Parameterisations of electron energy spectra and Cherenkov photon angular 
distributions in high-energy showers are developed, which serve as a basis for 
calculations of radiation emitted by air showers and may also be useful for 
other purposes. A direct application is the prediction of Cherenkov radiation 
as needed for reconstruction of longitudinal shower size profiles measured by 
experiments based on the fluorescence technique. 

The outline of the paper is as follows. In Sec. 2, an analytical expression is de- 
rived for calculating Cherenkov light production in EAS based on electrons 1 
only. In Sec. 3, electron energy spectra are shown to be universal, parame- 
terised in terms of a phenomenological definition of the shower age parameter 

1 Throughout this work, the notation 'electrons' include e + and e~ . 
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and applied for calculating the longitudinal Cherenkov light profile; the results 
are compared to other approaches and to a detailed CORSIKA simulation. 
Angular distributions of electrons and Cherenkov photons are investigated in 
Sec. 4. The former are shown to be universal and the latter are parameterised 
in dependence of altitude and shower age. Finally, the concept of ionisation 
energy deposit in simulation and reconstruction of fluorescence light profiles is 
discussed in Sec. 5. Part of the results presented here have already been shown 
in [12,18,21]. 



2 Calculation of Cherenkov light production 



The total number of Cherenkov photons diV 7 produced in a shower per interval 
of slant depth dX and angle d6 with respect to the shower axis is given by 



dA S(*,M) = E Ar&Ah) J y^KE) dE, 



dXd9 

^thr 



where dNi/dE (X) are the normalised differential energy spectra of corre- 
sponding charged particles of type % (= electrons, muons, ...) at depth X, 
El hr (h) the height-dependent Cherenkov energy thresholds, and y^h^E) the 
corresponding Cherenkov photon yields. A^X, 8, h) are the normalised angu- 
lar distributions of produced Cherenkov photons with respect to the shower 
axis. The photons are assumed to be produced at the shower axis, which is a 
good approximation for fluorescence technique applications since most charged 
particles are moving at a lateral distance of less than 60 m. 

The amount of so-called scattered Cherenkov light reaching a detector depends 
mostly on the total number of photons produced since the dominant Rayleigh 
scattering process is nearly isotropic. Therefore, the integral ^^-(X, h), corre- 
sponding to integration of Eq. (1) over all angles 9, is the most important quan- 
tity and typically used to calculate the scattered Cherenkov light contribution 
in fluorescence measurements. The so-called direct Cherenkov light contribu- 
tion, i.e. Cherenkov photons hitting directly the detector without scattering 
in the atmosphere, is determined by both the total Cherenkov light intensity 
and the angular distributions A % {X, 8, h). 

The number of Cherenkov photons produced by a charged particle of total 
energy E and charge Z in a wavelength interval between Ai and A 2 is given 
by 

dJVW, . 2itaZ 2 } 2 / 1 \ dA 
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where a is the fine-structure constant, (3 = v/c, p the air density at height 
h, and n the refractive index. Measuring dX along the shower trajectory we 
introduce the approximation that all particle trajectories are parallel to the 
shower axis, resulting in an under-estimation of produced Cherenkov photons 
of order one percent, cf. Sec. 3.4. The wavelength interval is usually about 
300-400 nm for fluorescence telescopes [1,5]. Since in air n ~ 1 and dispersion 
is negligible the integrand of (2) can be approximated as 

i - (M- 2 = i - (i - ^) _1 (i + s)- 2 - > (3) 



where 5(h) = n(h) — 1, and m the charged particle mass. The energy threshold 
condition for Cherenkov radiation in air reads 



E thT (h) = mc 2 /p5{h) , (4) 



giving 21 MeV for electrons and 4.4 GeV for muons at sea level. These thresh- 
olds increase to 37 MeV and 7.6 GeV respectively at an altitude of 10 km (US 
standard atmosphere). 

Particles heavier than electrons contribute less than 2% to the shower size 
around the shower maximum even for heavy primaries and their energy thresh- 
old for Cherenkov light production is very high. Therefore, to a good approx- 
imation, practically all charged particles can assumed to be electrons in the 
following considerations. Expression (1) simplifies to the following ansatz 

^l § (X,6,h) = A,(X,e,h)-N(X) J y,(h,E) f c (X,E)d\nE . (5) 

In E thl 



Here, N(X) is the charged particle number as function of depth X, and E thv 
the local Cherenkov energy threshold for electrons. For a given shower geom- 
etry, h = h(X) follows from the atmospheric model assumed, f e (X,E) is the 
normalised differential electron energy spectrum at depth X 

/ - M =5iijf)^ (x ' ¥) ' (6) 



and Ay(X, 9, h) is the normalised angular distribution of all Cherenkov pho- 
tons produced. 

We show in Sec. 3 and 5 that the ansatz (5) describes well the Monte Carlo 
results used as reference. 
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Fig. 1. Electron energy spectra obtained with CORSIKA for 3 different fixed 
shower ages. Shown are f 5 individual proton showers of 10 19 eV. 

2.1 Monte Carlo simulations 

The cosmic ray simulation code CORSIKA [19] version 6.137 with the hadronic 
interaction models QGS JET 01 [20] for high-energy interactions and GHEISHA 
2002 [22,23] for low-energy interactions has been used to study the develop- 
ment of high-energy showers. Energy spectra of electrons and angular dis- 
tributions of electrons and Cherenkov photons were obtained for individual 
showers and different combinations of primary energy, zenith angle and pri- 
mary particle type. The US standard atmosphere [24,25] has been used as 
atmospheric model. Also the longitudinal shower size profile and produced 
Cherenkov photon profile, which are needed to test the analytical model pro- 
posed in this paper, are calculated in the simulations. In CORSIKA calcu- 
lations, the wavelength dependence of the refractive index n(A) is neglected. 
This simplification is a good approximation as has been shown by a compar- 
ison [11] to GEANT3.21 [26] simulations taking into account the wavelength 
dependence. To reduce CPU-time, optimum thinning [27,28] of 10~ 6 has been 
applied unless otherwise noted. The UPWARD option has been enabled to 
follow the upward going electromagnetic particles 2 . Low-energy thresholds of 
100 MeV for hadrons and muons, and of 1 MeV for electrons and gamma-rays 



2 In CORSIKA, particles are counted at horizontal layers, which can be specified 
by the user. By default, particles with tracks at angles larger than 90° to the vertical 
are not followed. By applying the UPWARD option, these particles are treated in 
the electromagnetic shower component. 
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Fig. 2. Shower-to-shower fluctuations of electron energy spectra. The shaded bands 
indicate the RMS characterising the deviations of the individual showers of Fig. 1 
from the average. 

have been used, unless otherwise noted. 



3 Electron energy spectra 

In Figs. 1 and 2, normalised differential electron energy spectra f e (X, E) of 
individual showers obtained with CORSIKA are compared for fixed (effective) 
shower age 



' ' (1 + 2X max /X) ' v ; 

where X denotes the slant depth (in g/cm 2 ). The study covers the shower age 
range 0.8 < s < 1.2 that is most important for fluorescence observations. The 
spectra have been normalised according to 

1 dN r 
f e (X,E) = with J f e (X,E) dln£ = l, (8) 

where E cnt is the energy threshold adopted in the simulation (1 MeV in the 
examples shown). The importance of the low-energy threshold E cnt for nor- 
malisation and Cherenkov (and fluorescence) light calculations is discussed in 
more detail below. If the longitudinal development is described in terms of 
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Fig. 3. Universality of electron en- 
ergy spectra - different primary en- 
ergy. In shower age, mean showers of 
10 18 , 10 19 and 10 20 eV (each curve repre- 
sents the mean distributions of 15 indi- 
vidual showers) do not show significant 
differences in their electron energy spec- 
tra; here shown for s = 1.0 and proton 
initiated showers. 
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Fig. 4. Universality of electron en- 
ergy spectra - different primary par- 
ticle type. In shower age, mean pro- 
ton, iron and gamma-ray showers of 
energies > 10 17 eV (each curve repre- 
sents the mean distribution averaged 
over 10 18 , 10 19 , and 10 20 eV) do not show 
significant differences in their electron 
energy spectra, here shown for s = 1.0. 



shower age, the remaining shower-to-shower fluctuations are relatively small 
and due mainly to the depth of the first interactions. Deviations are small- 
est beyond the shower maximum because the individual distributions are the 
more insensitive to initial differences in shower development the older the 
shower is. The small fluctuations of the envelopes in Fig. 2 are of statistical 
kind and caused by the thinning applied in the simulations. The comparison 
is shown here only for proton showers of 10 19 eV for which the fluctuations are 
the largest and looks similar for different primary particle types at ultrahigh- 
energy (> 10 17 eV). 



3.1 Universality 



Average electron energy spectra f e (E,s) for proton showers of different pri- 
mary energies and showers of different primary particle types (proton, iron, 
and gamma-ray averaged over different energies 10 18 , 10 19 , and 10 20 eV) are 
shown 3 in Figs. 3 and 4. For the large electron energy range from 1 MeV up to 
10 GeV, the mean distributions are similar when normalised according to (8). 
The energy distribution of electrons in UHECR initiated EAS does not de- 
pend significantly either on primary energy or on primary particle type, which 
allows a universal parameterisation in shower age [16,18]. It is also indepen- 

3 In the case of gamma-ray showers, simulations were done for 10 19 ' 5 eV instead of 
10 20 eV as otherwise primary gamma-rays interact with the Earth magnetic field 
well before reaching the atmosphere (pre-shower effect) [29]. 
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E cut (MeV) 



Fig. 5. Dependence of shower size at 
maximum on the simulation energy 
threshold. The ratio of the mean shower 
size at shower maximum iV max to that 
extrapolated to E cut = is shown for 
different simulation thresholds. A pa- 
rameterisation of this ratio as a function 
of E cut is given [18]. 




10 '" 10 " 3 10 2 10 ' 1 1 



energy E [GeV] 

Fig. 6. Electron energy spectra obtained 
from CORSIKA (E cut = 50keV) with 
and without upward going particles. 
The spectra are normalised according to 
Eq. (8) assuming E cut = 1 MeV, for dis- 
cussion see text. 



dent of the shower zenith angle, however, this has only been investigated for 
inclinations smaller up to 60°. Increasing deviations from the spectral shape 
shown occur at 60° for high-energy electrons in the GeV-range [12]. 

3.2 Problem of low- energy cut and definition of shower size in simulations 

In general, simulating air showers is not possible without applying a low-energy 
cut -E cut on the tracked particles as the number of produced photons diverges 
for E cut — > [30]. The values of E cut that are typically applied in simulations 
are in the range of 100 keV - 3 MeV. Therefore, a parameterisation of electron 
energy spectra for the purpose of calculating Cherenkov light from simulated 
shower size profiles has to account for different energy cuts to be consistent 
with the simulation. Particularly in calculations of Cherenkov light based on 
ansatz (5), and for deriving energy deposit profiles to compute fluorescence 
light production, the low-energy cut applied in simulations is important since 
the number of charged particles provided by the simulations refers only to 
the particles above the threshold E cut . The effect of different values of the 
simulation threshold on the shower size is illustrated in Fig. 5. For a detailed 
discussion of energy deposit calculations, see [31,32]. 

According to the definition of shower size, particles are counted when crossing 
virtual planes (horizontal in the CORSIKA version used). Depending on the 
angular distribution, particles might not be counted when their tracks are 
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Fig. 7. Comparison of the new parameterisation Eq. (9), see also [18], and the energy 
spectra of an individual shower obtained with CORSIKA, proton, 10 19 eV. 

parallel to these planes and when they are going upwards. Multiple scattering 
of particles might cause low-energy particles (back-scattered) being counted 
multiply, which implies some ambiguity in the definition of shower size. This 
ambiguity is mostly related to low-energy particles as can be seen in Fig. 6. 
The simulated electron energy spectrum of a vertical shower is shown for 
both upward going particles having been accounted for and not. Differences 
in the corresponding distributions occur for electrons below about 1 MeV and 
increase with decreasing energy. Only electrons of less than about 1 MeV are 
affected because of their angular distribution (with respect to the shower axis), 
which is broad in the low-energy range up to a few MeV and steepens with 
increasing energy, see discussion in Sec. 4.1. 

3.3 Parameterisation in shower age 

Motivated by the high-energy limit of the energy behaviour of electrons in the 
cascade theory under approximation A [30] the following parameterisation 



is proposed. As shown in Fig. 7, the CORSIKA spectra can be reproduced 
well by ansatz (9) using the parameters given in the appendix. The parame- 
ter ao follows automatically from the normalisation condition (8). The ^cut- 
dependence of the normalisation ao is not negligible. In the energy range well 



fe (E, s) = a • 



E 



(9) 



(E + at)(E + a 2 ) s ' 



9 



10" 3 10" 2 10" 1 1 10 

energy E [GeV] 

Fig. 8. Mean electron energy spectrum derived with CORSIKA (circles) for s = 1.0. 
The Monte Carlo predictions are compared to the parameterisations according to 
Hillas [13] (dashed-dotted), the new parameterisation (Eq. 9) (solid) and Giller et 
al. [16] (dashed). 

below 250 keV, this dependence might be negligible, but for larger cut-off val- 
ues, the normalisation changes by up to e.g. about 10 % in the energy cut 
range from 50 keV to 2 MeV, see Fig. 5. A numerical expression for ao as a 
function of shower age s and E cut is given in the appendix. 

The parameterisation (9) is compared with other parameterisations and elec- 
tron energy spectra obtained with CORSIKA in the following. In the case 
of Hillas' parameterisation, the normalisation integral cannot be calculated 
directly. Given the fact that this parameterisation was obtained for 100 GeV 
primary photons, a larger disagreement to CORSIKA above 15 MeV (the lower 
validity limit given in [13]) might have been expected. Another functional form 
is proposed by Giller et al. in [16]. Neither Hillas' nor the parameterisation 
given by Giller et al. take into account different low-energy thresholds. While 
in [13] -E cut — 50keV is given, no energy threshold has been specified in [16]. 
In the following, the parameterisations of Giller et al. is normalised the same 
way as the parameterisation introduced here, see Eq. (8), using E cnt = 1 MeV. 
In the case of the integral distribution provided by Hillas, the parameterisa- 
tion is renormalised to compare properly. In this case the renormalisation is 
done according to the CORSIKA results on the dependence of shower particle 
content on simulation energy threshold, as displayed in Fig. 5. 

In Fig. 8 and 9, all three approaches are compared to the CORSIKA predic- 
tions. The comparison is shown for a mean shower averaged in shower age 
over 87 individual showers of different primary energies and particle types 
(10 18 , 10 19 , and 10 20 eV; 45 proton, 12 iron, and 30 gamma-ray induced show- 
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Fig. 9. Relative differences between model predictions and Monte Carlo results 
(relative to CORSIKA) as compared in Fig. 8 ([13] (dashed dotted), [16] (dashed), 
and the parameterisation (Eq. 9) (solid)). 

ers, cf. Fig. 3 and Fig. 4) at shower maximum. The other analytically calcu- 
lated distributions show a remarkable shift of the maxima with respect to the 
CORSIKA result of about 25 MeV at shower maximum (Hillas: ~ 40 MeV, 
Giller et al.: ~ 20 MeV), whereas the proposed parameterisation gives a very 
good overall description of the Monte Carlo results. The relative differences 
of the different approaches with respect to CORSIKA are shown for three 
stages of shower development in Fig. 9. The parameterisation (9) predicts the 
CORSIKA results within a few percent for the whole MeV-range, whereas 
uncertainties increase up to about 20 % at 10 GeV. In particular the energy 
range most important for Cherenkov light calculations (20 MeV to 1 GeV) is 
described with high accuracy. 

3.4 Total number of Cherenkov photons produced 

According to ansatz (5) the total number of produced Cherenkov photons 
can be calculated applying the parameterisation (9) of the electron energy 
spectrum f e (E, s) and the charged particle number N ch (X) provided by COR- 
SIKA. In Fig. 10, the result is compared to the CORSIKA simulated profile and 
to analytical calculations based on the other analytical formulae for electron 
spectra [13,16] already considered in the previous section. Using the analytical 
expression going back to Hillas leads to an over-estimation of the Monte Carlo 
result, whereas the work of giller et al. [16] results in an under-estimation. The 
calculation labelled 'Hillas (s fixed)' employs the parameterisation of [13] for a 
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Fig. 10. Comparison of the total number of Cherenkov photons produced per slant 
depth within an individual shower (proton, 10 19 eV, 30°) simulated with CORSIKA 
and analytically calculated by ansatz (5) using different parameterisations of elec- 
tron energy spectra. The different approaches are marked by dashed-dotted line 
for [13] and dotted line for the simplification s = 1 as it has been used by the 
HiRes/Fly's Eye collaboration [5,33], dashed lines for [16], and solid lines for the 
one proposed in this paper. The accuracy of Eq. (5) utilising parameterisation (9) 
is better than 2 % over the whole range important for fluorescence observations. 

fixed shower age of s — 1 only, as often used in shower reconstruction (see e.g. 
[5]). This approximation leads to a shift of the maximum of the Cherenkov 
profile by about 30-40 g/cm 2 towards larger depths, due mainly to the ne- 
glected reduction of high-energy electrons with growing age. Application of 
the parameterisation (9) reproduces the CORSIKA simulated number of pro- 
duced Cherenkov photons better than 1-2 % over the whole range important 
for fluorescence observations. 

The angular dispersion of charged particles, effectively increasing the Cherenkov 
yield per transversed depth dX along the axis, has not been taken into ac- 
count in the analytical approaches, cf. Sec. 2. Around the shower maximum, 
where the mean angle of electrons emitting Cherenkov light amounts about 
8°, the corrected yield would increase the predicted curves by about +1 %. 
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Fig. 11. Normalised electron angular distributions for an individual vertical proton 
shower of 10 19 eV. The electron angular distribution (with respect to shower axis) 
is shown for three different shower ages and various different ranges of electron 
energies. The electron angular distribution in high-energy showers is to a large 
extend independent of the shower age. 

4 Angular distribution of electrons and Cherenkov photons 



4-1 Electron angular distribution 



The energy dependent electron angular distribution with respect to the shower 
axis implies the angular distribution of resultant Cherenkov photons produced 
in EAS by convoluting in the Cherenkov emission angle. This angle slightly 
changes with altitude and amounts about 1.4° at sea level and decreases to 
about 0.8° at 10 km above sea level. Electrons undergo Coulomb scattering and 
the higher the particle energy the smaller is the scattering angle. Thus, the 
electron angular distribution is correlated with energy: The higher the mean 
electron energy, the smaller is the mean angle with respect to the shower axis. 

In Fig. 11, the normalised electron angular distribution with respect to the 
shower axis (averaged over azimuthal angle) is shown exemplarily for an in- 
dividual proton shower of 10 19 eV for three different shower ages, and various 
ranges of electron energies. The distributions do not change significantly with 
shower age within the statistical limitations. Deviations increase with increas- 
ing angles to the shower axis, especially for higher electron energies of a few 
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Fig. 12. Universality of (normalised) electron angular distributions. The electron 
angular distribution with respect to shower axis of numerous individual showers of 
different energies and primary particles is shown at shower maximum for various 
different ranges of electron energies, (proton, iron, 10 18 eV, 10 19 eV, 0°). 

hundred MeV and GeV. However, at that energies and angles, fluctuations are 
also caused by low statistics due to the thinning applied. That the electron 
angular distributions is practically independent of shower age supports the 
electrons energy to be the determining factor of their angular distribution. 

Charged particles are also deflected by the geomagnetic field. In the following 
we always average over the azimuthal angle and do not consider asymmetries 
implied by the geomagnetic field. Although such effect has been discussed 
not to make a gross change [13], the effect seems not to be negligible, an 
approximate treatment of the effect is given in [34]. At which extend the 
formulae given there describes the resulting azimuthal asymmetry sufficiently 
well would have to be investigated in detail separately. 

4-2 Universality 

As electron energy spectra in high-energy showers are universal and their scat- 
tering angle is mostly determined by the particle energy, the electron angular 
distribution should also be universal. 

In Fig. 12, the electron angular distribution of many individual proton and 
iron showers is shown for a fixed shower age (s = 1.0). For the MeV-range, 
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6 to shower axis [deg] 

Fig. 13. The angular distribution of Cherenkov photons with respect to the shower 
axis for 3 different heights obtained with CORSIKA is shown for one individual 
proton shower. 

the distributions of individual showers (of different primary energies and pri- 
mary particle types) do not differ much; larger fluctuations occur only for 
large angles. The angular distribution of electrons in the GeV-region shows 
a larger spread. This is in agreement with the universality studies of the en- 
ergy spectrum, where larger shower-to-shower fluctuations also occur in the 
GeV-range. 

In conclusion, the spectral shape of the angular distribution of electrons in 
high-energy showers is mostly determined by the electron energy and almost 
independent of the shower age (in the range considered of 0.8 < s < 1.2). It 
does not depend significantly either on primary energy or on primary particle 
type, see also [21,35]. 



4-3 Resultant Cherenkov photons 

The angular distribution A y (9,X,h) in Eq. (5) is the angular distribution of 
Cherenkov photons produced per angular bin with respect to the shower axis 
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normalised to one photon (averaged over azimuth): 

A,(6, X,h) ■ ^ = ^ (9, X, h) , with J A, d9 = 1 . (10) 

o 

Ay(0) follows from the angular distribution of underlying electrons. Thus, 
two dependencies follow for the angular distribution of produced Cherenkov 
photons in EAS: 

(i) Dependence on height 

The Cherenkov emission angle is slightly changing with altitude, this de- 
pendence is, however, negligible as the electron angles are much larger. Only 
electrons of energies above the Cherenkov energy threshold E thr are of in- 
terest, which implies that the higher E t ^ r the smaller is the mean angle 
of electrons emitting Cherenkov light. Finally, since the Cherenkov energy 
threshold E^ih) is - via the refractive index n{h) - a function of altitude, 
also the angular distribution of Cherenkov photons depends on height. This 
dependence on E thr can be approximated by an exponential function, see 
[34]. 

(ii) Dependence on shower age 

Since electron energy spectra develop with shower age, the part of electrons 
above the Cherenkov threshold energy E thr also changes with shower age. 
This is the reason why also the angular distribution of Cherenkov photons 
depends on shower age. 

Consequently, the photon angular distributions depend on both, the height 
due to the dependence of E^ih), and on the shower age because of the de- 
pendence of the electron energy spectrum f e {E, s). These dependencies of the 
angular distribution of produced Cherenkov photons can be seen in Fig. 13. 
With decreasing height and increasing index of refraction, the photon number 
at larger angles to the shower axis increases, whereas it decreases at angles 
smaller than about 2.5°, see also Fig. 14. 

4-4 Parameterisation in height and shower age 

In principle, one could use the universal electron angular distribution (finding 
a parameterisation in the mean electron energy or just use tabulated mean val- 
ues) for calculating the angular distribution of produced Cherenkov photons. 
As it turns out it is also possible to find a closed analytical approximation for 
the complete integral. 

It is common to describe the height dependence of the angular distribution of 
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Fig. 14. Angular distribution of produced Cherenkov photons with respect to the 
shower axis in a single CORSIKA shower for s = 0.8, 1.0, and 1.2. The Monte Carlo 
results are compared to predictions from [5] and the new parameterisation, Eq. (13). 

produced Cherenkov photons by an exponential function, where the scaling 
angle 9$ is a function of the Cherenkov energy threshold, see e.g. [5,34]: 

A 7 (h,9) = 1/0 O -e- 9/e ° . (11) 



Parameterisations of 9q as a function of local Cherenkov threshold energy 
Ethr(h) have been calculated by several authors; for example 9q = aE^, 
with (a, b) = (0.83,0.67) [36]. Traditionally this approximation is applied for 
calculating the Cherenkov contamination of fluorescence light signals from 
high-energy showers, see e.g. [5,33]. Generalising this ansatz to take into ac- 
count both effects, the dependence on refractive index as well as on the shower 
age we, write 

A^(9,h,s) = A y (s) ■ A y (9,h) 

= a s (s) ■ l/9 c (h) e-^W ? (12 ) 

where a s (s) is a polynomial of second order in shower age and the second, 
exponential term on the right depends on altitude only. The first accounts for 
so-called shower-to-shower fluctuations, namely the position of X max , and the 
latter includes the inhomogeneity of the medium, namely the refractive index 
changing with altitude. The most important range up to about 30° is described 
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Fig. 15. Quality of description of individual CORSIKA showers (iron of 10 18 , differ- 
ent zenith angles) by the new parameterisation Eq. (13). The statistical fluctuations 
due to thinning are reduced remarkably as very good thinning is applied (1CP 8 ). 
Showers are described practically independently of their inclinations. 

well by Eq. (12). To enlarge the range of validity up to 60° and improve the 
data description around 30°, ansatz (12) is extended by a second term to 



A,(6,h,s) = 0.(8) ^ e-W) + b a (s) ^ e- 9 ^ . (13) 

The numerical values of a s , b s , and 6 C , 6 CC , which have been found by a global 
fit to many individual showers of different primary particles and inclinations, 
are given in the appendix. The simulated Cherenkov photon angular distri- 
butions, which have been used to determine the coefficients of (13), although 
simulated with geomagnetic field, are averaged over azimuth. As shown in 
Fig. 14, the CORSIKA spectra are described properly by ansatz (13). The 
achieved quality of description by this ansatz is shown in Fig. 15. There, the 
relative differences with respect to individual CORSIKA showers of very high 
statistics (optimum thinning 10 -8 ) are shown. Taking into account both the 
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height and age dependence, also showers of different inclinations are described 
practically with the same accuracy of a few percent. 



5 Use of energy deposit for reconstruction and simulation of light 
profiles 



In fluorescence observations, one measures the fluorescence light profile, and 
the number of produced fluorescence photons in a shower is proportional to 
the local ionisation energy deposit dE^/dX: 

y (X) = y*(h) • ^(X) , (14) 



where is the fluorescence yield in air at altitude h. Thus, light profiles mea- 
sured using the fluorescence technique are measurements of the energy deposit 
by the shower rather than its shower size. The advantage of using dE dep /dX 
over shower size is not only that dE dep /dX is most closely connected to the 
measured fluorescence light but also solves the problem of the ambiguity in the 
definition of shower size in simulations as pointed out in Sec. 3. Furthermore, 
it can be conveniently simulated [31,32]. In addition, energy deposit profiles 
can be simulated with rather large low-energy thresholds E cnt > 1 MeV, and 
therefore is much less CPU consuming. 

For analytical description of the Cherenkov contribution to fluorescence pro- 
files, it is hence very useful to allow the calculation additionally as a function 
of d-Edep/dX, although the shower size is the physics quantity that is more 
closely related to Cherenkov production. 

5. 1 Parameterisation of mean ionisation loss rate 

We define the mean ionisation loss rate a e ff by 

a cS (X, E > E cut ) N ch (X, E > E cut ) = (X) , (15) 



where E cut is a low-energy threshold, which has to be applied in the case 
of shower simulations as discussed in Sec. 3. The ionisation loss rate is mostly 
determined by the charged particle energy. As a consequence of the universality 
of electron distributions, cf. Sec. 3 and 4, also a c s should neither depend on 
primary energy nor on particle type. 
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Fig. 16. The Mean ionisation loss rate a e fr as obtained from CORSIKA simulations 
according to Eq. (15) is shown for vertical showers of different primary energies and 
particles (10 17 , 10 18 , 10 19 eV, and proton, iron). The fit shown is given by Eq. (16). 

The corresponding study is shown in Fig. 16 where a e ff obtained from COR- 
SIKA simulations is shown for individual vertical proton and iron showers of 
different energies (10 17 , 10 18 , 10 19 ). In the simulations shown, E cut = 1 MeV 
has been applied. For the large range from about s = 0.5 to 1.2, the mean 
ionisation loss rate a e s(s) does not differ significantly for different primary 
particles and energies. The shower age dependence of a c s(s) can be described 
by 



a eS is) 



Cl 



(c 2 + s) 



+ c 4 + c 5 • s 



(16) 



which is the analytical expression superimposed in Fig 16; the parameter values 
are given in the appendix. 



5.2 Calculation of Cherenkov light as a function of energy deposit 



Using parameterisation (16) in shower age s, Cherenkov light production can 
be calculated as a function of energy deposit profiles following ansatz (5) 



diV 7 
dX 



(X,h) = N(X) J y 7 (h,E) / e (s,E) dln£ 



In E thT 
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Fig. 17. Comparison of total number of Cherenkov photons produced per slant 
depth within an individual shower (proton, 10 19 eV, 30°). Shown are the results 
of the CORSIKA simulation and analytically calculated by ansatz (5) and (17) 
respectively using the proposed parameterisation of electron energy distributions 
(Eq.(9)). " 



dE, 



dcp 



(X) J y 7 (h,E) f e {s,E) din E . (17) 



a eS {s) dX 

In E thr 

However, it should be noted that for proper calculations, f e (s,E > E cut ) and 
a e s(s,E > -Ecut) (and N(X,E > E cnt ) have to refer to the same low-energy 
threshold applied in the simulation. Moreover, vice versa, one can derive the 
energy deposit profile from a given shower size profile applying the parameter- 
isation (16). Such application is useful e.g. for simulating fluorescence profiles 
by so-called fast hybrid simulation codes like e.g. CONEX [37], which give 
shower size profiles as output. 

The Cherenkov calculation according to Eq. (17) is shown in Fig. 17 and com- 
pared to the CORSIKA result as well as to the calculation as a function of 
shower size. The same accuracy in reproducing the full Monte Carlo calcu- 
lation better than 2 % is achieved either using the shower size or the energy 
deposit profile provided by the simulations. 



6 Summary and conclusions 



For the purpose of developing an analytical description of Cherenkov light 
production in EAS, the universality of electron distributions in high-energy 
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showers has been investigated, namely the energy and angular distribution. 
In shower age, both have been shown to be independent of different primary 
energy and particle type to a good approximation. A parameterisation of the 
energy spectrum in shower age has been introduced, which describes showers 
independently of different primary energy, particle type and zenith angle at 
high accuracy of a few percent (within shower-to-shower fluctuations) for the 
electron energy range from 1 MeV to a few GeV, covering also the range most 
important for Cherenkov light emission. 

The electron angular distribution is mostly dominated by the particle energy 
and do not depend significantly on the shower age for a fixed energy. As a con- 
sequence, the angular distribution of produced Cherenkov photons has been 
parameterised as a function of height and shower age. The age dependence 
follows from the electron energy spectra. The dependence on height is due to 
the change of the Cherenkov threshold energy with altitude. 

Based on these universality features of high-energy showers, an analytical de- 
scription of the Cherenkov light production in EAS has been presented pro- 
viding both, the total number of produced Cherenkov photons as well as their 
angular distribution with respect to the shower axis. It offers the calculation 
of the direct and scattered Cherenkov contribution to measured fluorescence 
light profiles, see [38]. 

The advantage of using energy deposit rather than the number of charged par- 
ticles in the simulation and reconstruction of light profiles measured with the 
fluorescence technique has been pointed out. Having introduced a parameteri- 
sation of the mean energy deposit a e & (s), the scattered and direct Cherenkov 
light contributions to light profiles measured in fluorescence observations can 
be estimated either as a function of shower size or energy deposit profiles with 
an accuracy of a few percent. 

The geomagnetic effect leading to azimuthal asymmetries in electron and 
Cherenkov photon angular distributions seems not to be negligible and needs 
to be studied in more detail separately. 
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8 Appendix 



8. 1 Parameterisation of electron energy spectrum 



The parametrisation of the normalised electron energy spectrum in EAS is 
given by 

/c (E ' S) = ° ' (F + n WF-Un Vs ■ (18) 



The CORSIKA simulated electron energy spectra are described well by the 
following set of parameters: 



ai = 6.42522 - 1.53183 -s 

a 2 = 168.168 - 42.1368 • s , with E in MeV. (19) 

For applications in data analysis, the normalisation ao has been parametrised 
in shower age s for different E cut . The parameter ao is described by an expo- 
nential function in shower age 

a = k ■ exp(A;! • s + k 2 ■ s 2 ) , (20) 



where the parameters ko, k\ and k 2 calculated for six different threshold ener- 
gies in the typical range for simulation energy thresholds of 50 keV - 2 MeV 
may be linearly interpolated from the tabulated values given in Tab. 1 . 



Table 1 

Tabulated values for the normalisation ao (see Eq. (20)) of electron energy param- 





E cut [MeV] 


k 


h 


k 2 


2. 


1.48071e-01 


6.22334 


-5.89710e-01 


1. 


1.45098e-01 


6.20114 


-5.96851e-01 


0.5 


1.43458e-01 


6.18979 


-6.01298e-01 


0.25 


1.42589e-01 


6.18413 


-6.03838e-01 


0.1 


1.42049e-01 


6.18075 


-6.05484e-01 


0.05 


1.41866e-01 


6.17963 


-6.06055e-01 
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8.2 Parameterisation of angular distribution of Cherenkov photons 



The parameterisation of the angular distribution of produced Cherenkov pho- 
tons with respect to the shower axis is given by 



In this expression the age dependence is included by 



(21) 



a s (s) = a + ai • s + a 2 ■ s 2 (22) 
&s(s) = &o + ^i -s + b 2 -s 2 , (23) 

and the height dependence is taken into account by the expression 

6 c (h) = a ■ , with E thr in MeV (24) 

cc (h) = 7 • 6 c {h) , with 7 = a + ■ s . (25) 

The CORSIKA spectra are described properly using the following parameters 

(a , ai, a 2 ) = (4.2489 • 10"\ 5.8371 • 10 -1 , -8.2373 • 10" 2 ) 
(bo, h, b 2 ) = (5.5108 • 1(T 2 , -9.5587 • 10" 2 , 5.6952 • 1(T 2 ) 
(a, p) = (0.62694, 0.60590) 
(a, Pi) = (10.509, -4.9644) . (26) 

8. 3 Parameterisation of mean energy deposit 

The mean ionisation loss rate a e fr in shower age s can be approximated by 
cl 

a eS (s) = — + c 4 + c 5 • s , for E cnt = 1 MeV , (27) 

(c 2 + S) Cs 

with ci = 3.90883, c 2 = 1.05301, c 3 = 9.91717, c 4 = 2.41715, c 5 = 0.13180. 
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